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ABSTRACT 


A method of designing airfoils in cascade by means of 
conformal transformations is discussed. Design of these air- 
foils is regulated by five independent input parameters, with 
solutions obtained by digital computer. A large number of 
cascades generated by this method are compared. To evaluate 
the limits of performance, a parameter to indicate the ten- 
dency toward flow separation is introduced, with a limiting 
value established and verified. Proper solidity is shown to 
be of great importance in achieving low values of this separ- 
ation parameter. The value of proper solidity for a given 
blade thickness is shown to be relatively independent of 
fairing shape. To increase performance, reducing blade thick- 
ness with a corresponding increase in solidity is shown to 


yield much greater improvement than changes in fairing shape. 
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1.0 Introduction 

This study is an extension of a line of investigation 
originally conducted by Professor Theodore H. Gawain of the 
U. S. Naval Postgraduate School. The original investigation, 
sponsored by the Aerojet General Corporation, was specifically 
aimed at increasing the pressure rise per stage of axial flow 
pumps. The intended application was to reduce the size of 
propellant pumps for large liquid fuel rockets. The require- 
ments were for maximum pressure rise per stage, regardless 
of efficiency. This, in turn, demanded positive and accurate 
control of the pressure distribution around each blade, with 
particular emphasis on the adverse pressure gradient experi- 
enced by the boundary layer. 

Flow of a real fluid through an axial flow compressor or 
pump involves some formidable problems. Besides the effects 
of viscosity and turbulence; tip losses, centrifugal force 
fields, and other effects complicate the problem considerably. 
As is usually the case in engineering however, a simplifica- 
tion of the problem may yield useful results. If the stream 
surfaces of the flow are assumed to be concentric cylinders, 
then a cylindrical cross section of the flow may be repre- 
sented by an infinite linear cascade with two-dimensional flow. 
If all velocities are considered relative to the blade row, 
then the solution for the rotor and the stator have identical 
mathematical form. Following the usual aerodynamic practice, 


viscosity may be neglected outside the boundary layer. This 
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leaves only turbulence as a factor preventing the use of poten- 
tial flow considerations. In order to use the powerful analy- 
tic methods of potential flow, turbulence has been neglected. 
This does not mean however, that the data presented is invalid. 
Actually, the potential flow assumption is used to obtain the 
pressure distribution only. The relation between pressure 
distribution and separation was obtained from experimental 
data from actual turbomachines, and must therefore include the 
effects of turbulence. 

From the foregoing considerations, it may be seen that 
the study of cascade design has application to the whole field 
of axial flow turbomachinery. The prcblem of cascade design 
is much more difficult than the problem of isolated airfoil 
design. In the early days of airfoil design, much experimen- 
tal data was obtained on a bewildering variety of shapes. 
However, there was no reliable method of predicting the per- 
formance of an airfoil before it was tested. It was only 
after the NACA (now NASA) discovered how airfoil performance 
depended on shape, specifically on thickness, camber and fair- 
ing distribution, that any real progress was made in the field 
of airfoil design. 

The present state of the field of cascade design is that 
there is a large volume of experimental data available on a 
bewildering variety of cascades. Many efforts have been made 
to extrapolate the theory of isolated airfoils so that it will 


apply to cascades. However, at present there is no reliable 
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theory which covers the field of cascade performance in the 
way that NASA's theory covers the field of airfoil performance. 
The parameters of thickness, camber, and fairing distribution 
are still involved, plus a few more. The relations and in- 
teractions between these parameters are so complex, that it 

is very difficult to consider any aspect of cascade design 
separately. 

When considering cascades as complete entities, the inter 
action effects are all present, and need not be corrected for. 
However, experimental efforts in this direction have been nec- 
essarily limited in scope, due to the effort and expense in- 
volved in constructing and testing large numbers of cascades. 
Theoretical methods, such as the one outlined in this investi- 
gation, provide the means of considering complete cascades 
over a wide range of design parameters. This study was under- 
taken from the standpoint of discovering areas of possible 
improvement in compressor performance. 

Most of the efforts of the NASA in the field of compre- 
sor design have been for maximum efficiency with moderate 
output. In some cases, it may be required to attain maximum 
output, with perhaps moderate efficiency. The gains to be 
expected from these considerations are not of large order how- 
ever, aS a cascade develops its best efficiency not far from 
the point of stalling the blades. 

The basic problem, then, is composed of two parts. First 


it is necessary to devise a method for generating blade shapes 
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with the corresponding pressure distributions. For this pur- 
pose, the powerful methods of complex variables and conformal 
transformations are utilized. The method developed by Pro- 
fessor Gawain generates infinite cascades of airfoil shapes, 
based on five arbitrary parameters. Generally, the transfor- 
mations are similar to the Joukowski transformations with 
additional requirements. 

The second part of the problem is to establish a correla- 
tion between the pressure distribution around each blade and 
the occurence of separation. This problem must, of course, 
depend heavily on experimental data. Fortunately, work has 
been done in this area and some practical working limits have 
been established. See Reference 1. Based on this data, the 
simple parameter ( Vnax/V2)* has been chosen as an indication 
of the tendency toward separation. 

The author wishes to express his sincere appreciation to 
Professor Theodore H. Gawain of the U. S. Naval Postgraduate 
School for his assistance, guidance, and encouragement during 


the course of this investigation. 





2.0 The Method of Cascade Design 


2-1 Introduction 

Methods of designing airfoils by purely mathematical means 
have been known since the early days of aviation. Joukowski 
was the most notable pioneer in this field. He, and others, 
noted the possibility of designing cascades of airfoils by 
Similar methods. Since the early applications of cascade de- 
sign were to multi-winged airplanes, not much work was done in 
this field. As pointed out earlier, the application to turbo- 
machinery makes this field of considerable interest today. 

In considering the flow through cascades, the shape and 
position of any individual blade will affect the flow around 
that blade. In addition, the shape and position of all the 
other blades will also affect the flow around that same blade 
at the same time. rhea hacen more parameters must be used in 
cascade design than are used in isolated airfoil design. In 
addition to blade shape, the relation of the blades to each 
other, and their relation to the flow must all be described. 
The method developed by Professor Gawain uses five independent 
parameters to describe an individual cascade. 

The solution for blade shape and velocity distribution 
corresponding to chosen values of the five parameters is found 
by using relations of potential flow and conformal mapping in 
three separate planes. These have been termed the near circle 


plane, the circle plane, and the cascade plane respectively. 


a, 
__ 
re ee ee 
(=~ © =f ae oe 1a 
Ae Fn me eee oe —— ~— 
——t EO a omen voles 
= é i fain 


“= & 





In reading through the following sections, which describe 
the general method of generating cascades, reference may be 
made to Figure 5, which illustrates a typical contour in each 
of the planes referred to above. A summary of the mathematical 
development of the method is given in Appendix A. A more de- 
tailed analysis may be found in Reference 2. 

The solution of these equations, many of which are trans- 
cendental, is an impossible task by hand calculation. The 
entire problem has been arranged in FORTRAN language for solu- 
tion by a digital computer. The CDC 1604 computer at the U. S. 
Naval Postgraduate School was used for this investigation. 

The basic computer program, along with later modifications is 
included in Appendix B. 

The following sections describe the general method of 
generating cascades from five input parameters. Each parameter 
is described in detail as it enters into the calculations. 


The parameters are: 


/8 = Mean flow direction 

6, = Turning parameter 

Oo” = Equivalent flat plate solidity 
& = Thickness parameter 

K = Shape or symmetry parameter 


2.2 The Flow Parameters 
Flow through a cascade is usually described in terms of 


two velocities. The first velocity is that far enough up- 
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stream of the cascade so that local effects due to the cascade 
are not evident. This is called the inlet velocity, and the 
angle it makes with the reference axis is called @,. Similarly, 
the velocity far enough downstream so that local effects have 
damped out is called the outlet aap eye and its angle with 
the reference axis is called #3. The magnitudes of these two 
velocities are not independent, since the laws of continuity 
require that the component perpendicular to the cascade axis 
be equal for the two cases. This concept allows velocities to 
be conveniently nondimensionalized, so that only angles are 
important. 

Figure 1 is an illustration of two typical cascades. An 
appropriate velocity vector diagram is shown for each. Note 
that only two parameters are necessary to describe the complete 
diagram. Instead of the angles G) and G2 it was more conven- 
ient to use two different parameters. Consider the following 
expressions: 

Tan @ = 4(Tan A, + Tan G2) 
B, = (Tan /, -Tan £2) 
The first expression defines the angle (, or mean flow direc- 
tion, which is the first input parameter. The second expres- 
sion defines the turning parameter B,, which is the second 
input parameter. The relationship between all these quantities 
is shown in the diagram on Figure l. 


The starting point of the calculations is the near circle 
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plane. In this plane, the complex potential for a series of 
doublets spaced at intervals of 2 /J along the y axis is conm- 
bined with the potential for uniform flow at the angle / with 
the x axis. The velocity of this flow is set at unity, which 
then becomes the reference for all the nondimensional veloci- 
ties derived later. If only one doublet were used, the resul- 
ting potential would describe the flow without circulation 
around a single circular cylinder. In this case, with an in- 
finite series of doublets along the y axis, the resulting 
potential describes the flow without circulation around an 
infinite series of nearly circular bodies spaced along the y 
axis. These bodies are not actually circles because the pre- 
sence of the infinite stack of doublets distorts the flow 
around any individual doublet, such that the flow contour is 
a slightly flattened circle. The amount of distortion is re- 
lated to the strength of the doublets and the spacing between 
them. In this case, the spacing between doublets is fixed at 
277 , so doublet strength is what determines the relation 


between spacing and size of the bodies in the near circle plane. 


2.3 The Solidity Parameter 
The strength of the doublets is fixed by the input para- 


meter oO-~ , or the equivalent flat plate solidity. This para- 
meter is closely related to, but not identical with, the ordinary 
solidity c/S, or the chord-pitch ratio of the blade row. The 


distinction can best be illustrated by an example. Generally, 
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any row of arbitrary blades of solidity c/S may be transformed 
conformally into a related row of flat plate vanes. The chord- 
pitch ratio of the latter is then denoted by @ . It follows 
that if two arbitrary cascades have the same value of © , 

then either may be conformally transformed into the other. But 
if the value of © is different for each, then no such trans- 
formation is possible. © then, may be regarded as a kind 

of generalized solidity, of perhaps deeper aerodynamic signi- 
ficance than simple geometric solidity. 

The relation between values of © and geometric solidity 
is directly influenced by blade thickness. Very thin blades 
have geometric solidities almost exactly equal to © . With 
thicker blades, the geometric solidity is always less. The 
reason for this may be illustrated by examples of two extremes. 
If a contour in the near circle plane is transformed confor- 
mally into a flat plate, the chord length of the flat plate 
will be nearly twice the diameter of the original contour. 
However, if no transformation is made at all, an "airfoil" of 
100% thickness results, whose chord length is equal to the 
diameter. Both cascades would have the same value of O~, but 
the geometric solidity of the latter case would be half that 
of the former. Between these two extremes, finite thickness 


airfoils in cascade have solidities less than G7 . 


2-4 The Turning Parameter 


The turning parameter, B,, is a measure of the change in 
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direction of the flow passing through the blade row. It is 
actually a nondimensional form of the downwash velocity induced 
by the blades. There is a very simple relation between B, and 
the circulation around each blade. In fact, for a constant 
mean flow angle, this relation is a simple direct proportion. 
This convenient state of affairs is the reason why the two 
flow parameters were expressed in their present form. 

In isolated airfoil theory, the method of introducing lift 
or circulation on an airfoil, is to add to the complex poten- 
tial, a term representing a vortex of suitable strength. How- 
ever, mere addition of a vortex will produce a distortion in 
the contours in the near circle plane, since these contours 
are not exactly circles. In addition, an infinite stack of 
vortices is necessary to produce circulation around each of the 
blades in the cascade. To avoid these resulting undesirable 
perturbations in what will ultimately become the airfoil shape, 
and indirect method of adding circulation is used. 

In this method, a transformation is made from the near 
circle plane to the circle plane. An equation is developed 
which transtorms the flow in the near circle plane into the 
flow around a perfect circle, centered on the origin. The cir- 
culation term may then be introduced without fear of distor- 
tion. The transformation is then made back to the near circle 
plane. Since the same equation is used to transform both ways, 
the original flow contour remains unchanged; however, the 


velocities at each point are now different, due to the circu- 
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lation. 


2.5 The Blade Shape Parameters 


The final transformation is trom the near circle plane 
to the cascade plane. This is done by the Joukowski method 
of shifting the axes slightly in such a way that the flow con- 
tour transforms into an airfoil shape. Two arbitrary para- 
meters were chosen to control the blade shape. Cne of these 
is called the thickness parameter, and is given the symbol é€ . 
This quantity is generally related to the thickness, and also 
to nose curvature, but not in any fixed relationship. Large 
values of € produce thick airfoils, and small values produce 
thin airfoils. Figure 2 is an illustration of several airfoils, 
showing the effect of & . The limiting value of O for € does 
not however produce an airfoil of zero thickness. Due to the 
Slight distortion of the contours in the near circle plane, 

&=0 may produce the mathematically significant, but physi- 

cally useless case where the "lower" surface of the airfoil 
crosses through the "upper" surface. € then, can be small but 
not zero, and the airfoils produced can be thin but not of zero 
thickness. 

The other parameter is the shape or symmetry parameter K. 
The effect of variations in K may be seen in Figure 3. Gener- 
ally, K=O produces a blunt leading edge, and cusped trailing 
edge similar to a Joukowski airfoil. K=l will produce a sym- 


metrical shape which is generally a cambered ellipse. As K 
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increases from zero to one, the point of maximum thickness 
moves from the quarter chord to the half chord point, and the 
trailing edge rounds off. Values of K larger than one will 
produce airfoils facing the wrong way (Trailing edge into the 


flow) and hence must not be used. 


2.6 Stagnation Points 


During the transformations made in this analysis, the 
stagnation points serve as reference points to locate the re- 
lation of the contours in the various planes. For this reason 
there is always a stagnation point at the leading and the 
trailing edge. Having a stagnation point at the leading edge 
is sometimes called the Theodorsen condition. This means that 
in all the cascades developed by this method, the airfoils are 
operating at the ideal angle of attack. Variations in lift 
then, are achieved by variations in the camber of the blades. 
Figure 1 graphically illustrates this fact. Two typical 
cascades are drawn, with the principal variation being in turn- 


ing parameter, B The values of C were adjusted to give 


6° 
the same blade thickness. The difference in camber may be 
clearly seen. 

The condition of a stagnation point at the trailing edge 
is called the Kutta condition. This fact has particular sig- 
nificance when considering the blade shapes where K=l. As 


pointed out previously, these shapes are generally cambered 


ellipses. It is a well known fact that flow around an ellipse 
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of a given thickness ratio produces a lower velocity peak than 
the flow around any other shape of the same thickness ratio. 
This class of blades then should exhibit the lowest loading 
attainable for any possible fairing shape. Of course, any 
attempt to duplicate this fairing shapein an actual cascade 
is doomed to failure because of the rounded trailing edge. 
However, this case still has significance in two respects. 
First, the performance of a cambered ellipse, with the Kutta 
condition mathematically imposed, represents a useful limiting 
cease of performance due to fairing shape. Secondly, a blade 
Similar to this, with modifications to produce a sharp trail- 
ing edge, would have very nearly the same performance. The 
critical area of flow near the leading edge would not be af- 
fected much by modifications at the rear of the airfoil. 
Figure 4 is a plot of the pressure distribution around 
the airfoils of Figure 3. The above premise is verified, as 
the negative pressure peak can be seen to decrease as the 
value of K is increased. Note that the pressure distribution 
aft of the peak is essentially linear to the trailing edge. 
Reference 3, p- 354 states that this type of pressure distri- 
bution is typical of decelerating cascades. Note also that 
the slope of this line becomes less as K is increased, indi- 


cating a less adverse pressure gradient in the boundary layer. 
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3.0 Prediction of Separation 


3.1 # Pressure Gradient 

The phenomenon which ultimately limits the attainable 
pressure rise of a cascade is separation. Flow separation 
occurs when flow in a boundary layer experiences an adverse 
(positive) pressure gradient. This adverse gradient tends to 
decrease the velocities in the boundary layer. Ultimately, 
the point is reached where the velocities close to the surface 
have been reduced to zero. At this point, flow separation 
occurs. Beyond this point, some velocities become negative, 
and a region of reverse flow exists. The resultant effects 
cause high drag and loss of circulation around each blade 
element. Many experiments have verified that any adverse 
pressure gradient, however small, will cause separation if it 
acts long enough. The concept of a gradient and the distance 
over which it acts, leads one naturally to suspect pressure 
ditrerence as a significant parameter in predicting separation. 
The gradient, integrated over the distance is mathematically 
equal to the difference in pressure. It was previously shown, 
furthermore, that the pressure distribution over typical air- 
foils in this system is approximately linear aft of the nega- 


tive pressure peak. 


3.2 Velocity Considerations 


From another point of view, one might examine the velo- 


cities of flow through a cascade. From the inlet velocity 
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described previously, velocity is prought to zero at the 

front stagnation point, accelerated to a local maximum some- 
where on the surface of the blades, reduced toward zero at 

the trailing edge, and finally settles down to the outlet 
velocity far behind the cascade. Besides the inlet and outlet 
velocities, only one other velocity is really significant, and 


that is the local maximum, V From these three quantities, 


max ° 
only two independent dimensionless parameters may be construc- 
ted. If the inlet velocity and Vmax Were used, this parameter 
would be related to the flow around the forward part of the 
airfoil. Since this flow involves a negative pressure grad- 
ient, it is probably not very significant with respect to 
separation. Similarly, the quantities Veex and the outlet 
velocity, called Vo» are related to the flow around the rear 
part of the airfoil, which involves a positive pressure grad- 
jent. These two parameters then, combined in some dimension- 


less manner, should act as an indicator for the tendency 


toward separation. 


3.3 Diffusion Factors 
From its many experimental studies of flow through cas- 
cades, the NASA has arrived at essentially the same conclusion. 
They have used several dimensionless combinations of — and 
Vo to define blade loading limits. One of these is the 


"Local Diffusion Factor" which is: 
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max 
This factor is described in Reference 1 and supporting data 
indicates that separation is likely to occur for values of Dy 
greater than 0.5. 

Since the velocity distribution around the blades must 
be known to find values of Dy , it is necessary to use experi- 
mental data, or an analytical treatment such as the one des- 
cribed in this investigation. To get around this difficulty, 
the NASA also defines an approximation, D, which may be ob-= 
tained from more easily determined quantities. D, is called 


the "Diffusion Factor" and is given by: 


D= 1- CoS A) + COS Ai (Tan 71 - Tan £2) 
Cos 785 2% 
It is seen that this quantity can be determined from the inlet 
and outlet velocities, and the solidity of the blade row. The 
supporting data indicates that separation is likely to occur 
for values of D greater than 0.6. 

Since the methods used in this investigation furnish all 
of the flow quantities described above, as well as the com- 
plete velocity distribution around each blade, the opportunity 
to compare the two NASA parameters is presented. As will be 
described later, the correlation between these two parameters 


is not too good for the solidities considered here. 
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3-4 The Boundary Layer Loading Parameter 


From the above development. it may be seen that the para- 
meter to indicate separation should involve the pressure ratio 
across the rear part of the airfoil, and be some dimensionless 
combination of the quantities Vo and Viiax° The NASA has in- 
dicated that most compressors operate in the range where 
Reynolds number effects are insignificant. Therefore, it is 
felt that no particular inaccuracies result from omitting 
consideration of Reynolds number. 

In order to avoid using one parameter for pressure dis- 
tribution, and another for the tendency toward separation, 


the quantity S was selected where: 
= 2 
- ( Vinax/ Vo) 


This quantity is referred to as the boundary layer loading 
parameter, and is quite versatile in application. It is 
directly related to the pressure coefficient, and in fact only 
differs from that quantity by unity. S also has a constant 
functional relationship with the NASA parameter D). The value 
S = 4.25 corresponds to the value of D, of 0.5, and has been 
selected as the limit loading in the boundary layer for blades 


in this investigation. 
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4.0 Results and Discussion 

Theory and experiment both point toward increasing 
solidity as the method of obtaining increased pressure rise 
per stage. As solidity is increased, the lift coefficient on 
each blade is decreased, and hence the loading on the boundary 
layer is decreased also. However, as the solidity is increased, 
there is an adverse effect due to the blockage caused by the 
finite thickness of the blades plus the boundary layer. 
Finally, the point is reached where the beneficial effects of 
decreased blade loading are just counterbalanced by the ad- 
verse effects of blockage. The quantitive determination of 
this point was one of the major objectives of this investiga- 


tion. 


4.1 Organization and Scope 


Since actual cascades are designed on the basis of blade 
thickness and geometric solidity, it was felt that the results 
of this investigation ought to be organized on this basis, 
rather than on the basis of the mathematical parameters & and 
Oo” . The value of blade thickness is not directly predictable 
from the input parameter € . Therefore, the basic computer 
program was modified to use an iterative scheme to achieve the 
desired values of blade thickness. Details are given in 
Appendix B, which describes the basic computer program and the 
modifications. Seven blade thicknesses were selected to cover 


the field of possible blades. These ranged from the limiting 


18 






















—_— Pe 
eC — ee  / 
- 
2°" > = ee os | 
ore 4.——————— ea @ tenes 
eer a 


ee + = et | 








———_-— ll 


- 8” “=~ =m hii») iii ZZzz-, 
See at —— ry Ws 


SL ————— —_ CC ae 


——- « — ——— 


—_—_—-- -<s — jp 

6 & io ——— —»-s eaten * 
eeatlien a ee, 24 
—! = ‘ -_ —_—_ mle 
° — — we - eo, wl. 


i 
« 
- 
ao 


case of extremely thin blades to ~<4% thickness blades in four 
percent increments. 

It was decided to conduct this investigation using only 
one value of average flow angle . This value was chosen to 
be 45°, which was a somewhat arbitrary decision, but not with- 
out some sound theoretical basis. It is stated in Reference 
4 that maximum efficiency may be obtained when the average 
flow angle is near 45°, and many compressors are built accord- 
ingly. 

In spite of reducing the parameters to four by the above 
decision, over 1470 separate cascades were generated and ana- 
lyzed for this investigation. Additional parameters could 
have been introduced to give more control over the blade shape, 
but they would have greatly magnified the problem. As the 
data will show, additional blade shape parameters would not 


have affected the significance of the results very much. 


4.2 Series 1 Plots 

The first plots of data from the computer run were made 
ef boundary layer loading vs. solidity, for airfoils of con- 
stant blade thickness. Five values of the shape parameter K 
were used, and data was obtained at six values of turning 
parameter, B,. This meant that there were 30 graphs in series 
one. Figures 6 through 9 are representative graphs from this 
series. 


It was noted with satisfaction that the curves showed a 
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pronounced minimum loading as soliaity was increased. The one 
exception was the case of the extremely thin airfoils, for 
which the curve appeared to be asymptotic to some low value of 
loading. This is not unreasonable, since the blockage effect 
would be very small for this case. The case of thin airfoils 
then, representSa sort of limiting value for cascades. The 
point at which minimum loading occurred for a given thickness 
airfoil was called the optimum solidity. 

It may be noted that while the value for loading at an 
optimum point could be determined with accuracy, the corres- 
ponding solidity for that point was subject to some uncertainty. 
Therefore, another point was introduced where the solidity 
could be determined with greater certainty. This was called 
the design point, and is defined as the point toward lower 
solidity where the loading equals 1.10 times the minimum value. 
There is further justification in this step from the fact that 
the basic analysis neglects boundary layer. The blockage 
effect of a blade plus boundary layer must be greater than the 
effect due to the blade alone. Therefore, the true optimum 
solidity must be somewhat less than the apparent optimum based 
on non-viscous flow. The exact relation cannot be found with- 
out a detailed analysis of the boundary layer. This simple 
method sidesteps this formidable problem, and is at least 
qualitatively in the right direction. 

For the case of higher values of turning parameter and 


thick blades, there was some scatter in the data. Since the 
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data is all analytical, and should therefore be reasonably 
regular, some effort was made to find the source of the scat- 
ter. The relationship between the parameters © and ¢/S5 was 
suspected first. Comparative plots of these two parameters 
however, showed that they varied in a regular fashion with 
blade thickness. It was then decided that the scatter was due 
to the fact that velocities were measured only at 40 points 
along the airfoil. If the true maximum velocity occurred be- 
tween these finite points, then the value for boundary layer 
loading parameter would be in error. Remembering the fact 
that this error could be only on the low side was an aid in 
fairing the curves through the scattered points. 

It was noted that the data became more regular as the 
value of K approached one. For these airfoils, either the 
flow was more regular about the nose, with no sharp peaks in 
velocity, or the point of maximum velocity happened to coin- 
cide more often with one of the finite points. Examination 
of the pressure plots of Figure 4 points toward the former 
as the cause. With high values of K, the pressure changes are 
seen to be more gradual, and of lower magnitude. Because of 
this, the actual maximum velocity was always close to the 
velocity recorded at one of the 40 points. 

Figure 10 is an illustration of the effect of K in reduc- 
ing the boundary layer loading for a given value of turning 
parameter. With all other factors equal, the loading steadily 


decreased as the value of K increased toward unity. Note that 
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the effect of K was much greater tor the thick blades. For 
the case of extremely thin blades, K had no effect at all, 


which was not an unexpected result. 


4.3 Series 2 Plots 

Since the original series involved 30 graphs, some effort 
was made to summarize this data. The result of this effort 
was Figures 11 through 15, which collect all the data for a 
particular value of K on a single page. Boundary layer load- 
ing parameter and solidity are plotted along the coordinates 
as before, but now the solidity is design solidity. In effect, 
all the design points for a particular value of K have been 
plotted. Points of constant values of turning parameter have 
been connected, as have points of constant airfoil thickness. 
The result is that four interdependent parameters are shown 
for each point on the chart. Any chosen point then, repre- 
sents a particular value of blade thickness and turning para- 
meter, with the corresponding boundary layer loading and 
design solidity for the cascade. 

Figures 16 through 20 are plots made in exactly the same 
manner, except that optimum points were used instead of design 
points. It was necessary to change the scale for solidity 
along the horizontal axis, since optimum points represent 
higher values of solidity in every case. In addition some 
difficulty was experienced with irregularity in the data. 


This was due at least in part to the uncertainty in determin- 
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ing the value of solidity at the optimum point. This optimum 
point data was included as a limiting case, to illustrate the 
maximum solidities to be used when the boundary layer is neg-~ 
ligibly thin. 

The choice of solidity as a coordinate is slightly mis- 
leading. It might be inferred from this choice that solidity 
should be used as a design input. However, solidity was 
plotted in this manner purely as a convenience in transferring 
data from the series one graphs. In design work, the two 
"inputs" could be boundary layer loading, and turning para- 
meter. The “output” information would then be blade thickness 
and solidity. An attempt was made to plot boundary layer 
loading parameter vs. turning parameter, but the resulting 
blade thickness and solidity lines were so nearly parallel 


that the information was difficult to read. 


4.4 Comparison of NASA Separation Parameters 


The data on the series 2 plots contains all of the infor- 
mation necessary to compute both NASA separation parameters. 
This affords an opportunity to compare these two parameters 
directly. The boundary layer loading parameter S, and the 


NASA parameter Dy have a constant relationship such that: 


Thus a Dy equal to 0.5 corresponds to a value of S of 4.25. 


This value is entered as a heavy line on all series 2 plots. 


23 










i aa 
eS ee 
Sy 2 : 












a 
» 
m- 








Storia 
a ee 
ee 
7 eee! 6 mare + A Se hme 2 
Se aery or suse 
ss ss ——- 
—> —)* + ee o 
et eee te = ee 
Celene - | Se. : 

: Pee tw 6 coe bees pe Ae 
















ss 








= = 
ad 

































—~*> = 
=? 4 
<= — ly — 
———— tte OGD 66 ee 848 Shee | 
ee —_—— na 


> = —— =) — = ec 





21—& es | —a + inte?) ot 


—~ 6 —TE ee S & -« — i + a | 





The approximation LD, is relatea to the parameters B, Bo» 


and c/S by the following relationship: 


B, 


a 


ELS 1 4 (tan/3+ +) 


Thus, since (2 has been fixed at one value, each value of By 





and c/S will result in a particular value for D. This expres- 
sion for D shows that it cannot be a good approximation over 

a wide range of solidities, since it indicates that loading 
will always decrease whenever solidity is increased. It has 
already been demonstrated that as solidity increases beyond 
the optimum point, loading increases also due to the blockage 
effect. 

Lines corresponding to a value of D=.5 and D=.6 have been 
plotted on Figures ll through 20 as dashed lines. A note of 
caution is necessary in using the data on these figures. It 
must be remembered that there is really only one solidity 
shown on each chart, and that is the design (or optimum) soli- 
dity. The fact that solidity is shown as a variable is due 
solely to the fact that each combination of blade thickness 
and turning parameter demands a particular value of solidity 
at the design (or optimum) point. 

If the parameter D is a true approximation for D), then 
the lines corresponding to D=.5 and Dj=.5 should be nearly 


the same. Examination of the figures shows that there is a 


2k 






<venot 
Sere (i 
“4 + ari! 4 













a — _ 





=> 

















ae pe eet at 
ene “pata 
— 2 ey oe ee 
Care) tee) eee) nt ee ott ie & 
i ts pe le cence eel 
ed ree Fe ee ee! eee ee 
tmnt! eh bee Ay 
























Sa 
j 7 t-« _ 
> @ 2 
























ee et a ee ee de et eye ee 
=p + ml! eee a eee” 4 pen 
0 tte el el es 1) eee 
ee te Ale Bem me grey tem —ee fs 
. Ne i  L 
2 ieee ere 6 ee ee ot) 
SA el be es ee cee! ned oe pe Hole 
wml > ane’ — | 
imi ar tie a = 

ee re rryee ew me) mo ae al 
ae es tie fh ee) ee eee ere auneg ala 
7 eT eee al a! 





















fair match for the case of optimum solidity. Generally, D 
gives a value of loading which is nigher than D,, but the 
error is only around 10%, which is not unreasonable. For the 
case of design solidity, there is a greater discrepancy. For 
very thick blades, D still imaiicatos a value of loading higher 
than D), but this relation reverses so that for thin blades 
Dy indicates the higher loading. 

The data in Reference 1 states that separation is likely 
to occur for values of D) greater than .5 or for values of D 
greater than .6. This in itself is an indication of a mis- 
match between D and Dj). If the occurrence of separation is 
taken as the matching criterion, then the line corresponding 
to D=.6 should match the D)=-5 line. For the case of optimum 
solidity it may be seen that the matching is very poor, where 
D indicates a much higher loading than D,° At design solidity 
the matching is better, but D still indicates values of load-~ 
ing which are around 18% higher than D,). Extrapolating this 
trend leads one to suspect that, for some value of solidity 
less than design solidity, the two parameters D=.6 and Dj=.5 
would probably match quite closely. From this, it is not un- 
reasonable to expect that the NASA obtained its experimental 
data from cascades whose solidity was less than the design or 
optimum as defined in this report. A few checks of source 
data for Reference 1 have verified this premise. Solidities 
were generally 1.0 for cascades with blade thicknesses of 


eight and ten percent. A check of the series 2 plots will 
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reveal that this is less than the design or optimum point for 


any fairing shape. 


4.5 Comparison of Blade Thickness and Fairing Shape Effects 
Figures 11 through 20 clearly show the gains that can be 


made in turning parameter, at no increase in loading, by 
using thinner blades at the corresponding higher solidities. 
However, it was desired to compare the effects of changes in 
fairing shape and blade thickness on the same chart. Figure 
21 was made using the design point data, and Figure 22 was 
made with data from the optimum points. The value of S = 4.25 
was selected as the best estimate of maximum permissible blade 
loading. The values plotted were turning parameter vs. solid- 
ity at this constant value of blade loading. This resulted in 
lines which represented varying blade thickness at a constant 
value of K. Lines of constant blade thickness were then drawn 
in. It is interesting to note that these lines are nearly 
vertical. This means that the value of design or optimum 
solidity for a particular blade thickness is relatively inde- 
pendent of blade shape. This result was not directly predic- 
table from theory, and may prove useful in planning experi- 
mental investigations. 

As blade thickness is decreased, it is clearly seen that 
the importance of fairing shape decreases correspondingly. 
When blade thickness reaches about four percent, there is 


almost no discernible difference due to fairing shape. To 
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compare the relative effects of tuiring shape and thickness 
changes, consider the following example. Select a blade thick- 
ness of 12% with K = .010, which is a relatively poor fairing 
shape with a blunt nose. For this shape at the proper solid- 
ity, a turning parameter of .74 is possible without separation. 
If K = .990 the shape is nearly the best possible. 

Making this change increases the permissible turning 
parameter to .83, or a 12% improvement. This same improvement 
could be obtained with the original shape by decreasing the 
thickness to a little less than 8% and increasing the solidity 
slightly. If the blades could be decreased to 4%, with the 
corresponding increase in solidity, the turning parameter 
could be increased to .92 without changing the loading. This 
is a 24% improvement, or double the best improvement that 
could be obtained by fairing shape changes alone. As a limit- 
ing value, if the blades could be made infinitesmally thin, 
the increase in turning parameter could be as high as 50%. 

It must be remembered that these performance figures are 
for cascades operating only at the ideal angle of attack. It 
may well be that these increases in peak performance are at- 
tained with some sacrifices in flexibility. This would be of 
no consequence in a constant output device, but for general 
application, the problem of operating off-design must be con- 
sidered. The analytical method described in this study can 
be modified to consider the case of operating at angles of at- 


tack other than the ideal. This matter is discussed more 


fully in the section on suggestions for further work. 
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5.0 Conclusions 

It is concluded that the method of generating cascades 
as described in this report is a workable one, and provides 
useful results. Cascades are generated as complete entities, 
which eliminates the need for separating the effects of the 
various design parameters. Use of a high speed digital com- 
puter means that large numbers of cascades can be generated 
and analyzed quickly. This allows exploration over a wide 
range of design parameters. 

solidity in experimental cascades has often been arbi- 
trary. This investigation has revealed that the matter of 
solidity is of paramount importance. For any combination of 
fairing shape, thickness and camber, there is a particular 
value of solidity which will result in the lowest loading of 
the boundary layer. This loading can increaséquite rapidly 
as the solidity is varied from the proper value. 

Data obtained in this investigation has yielded the un- 
expected result that the optimum value of solidity referred 
to above is relatively independent of the fairing shape of the 
blades. An important application of this result will be in 
planning experimental investigations of cascade performance. 
The actual values of solidity given in this report may be used, 
with the probability that they will be close to the optimum 
values for any sort of physical blade shape. 

Three parameters intended as indicators of the tendency 


toward flow separation were compared in this report. Asa 
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result of this comparison, it is concluded that the parameter 
S, as defined in this report, is a useful indicator in this 
respect. It is further concluded that the value S = 4.25 is 
an acceptable design limit, which is as consistent with NASA 
data as this data is consistent with itself. 

With respect to the problem of increasing the performance 
of airfoils in cascades, it is concluded that the matter of 
blade thickness is of much greater importance than the fairing 
shape. The example was given of a cascade with blades of 12% 
thickness. Possible increases in turning parameter for this 
cascade due to improvements in fairing shape were shown to be 
on the order of 12%. Increases in turning parameter produced 
by decreasing the blade thickness to 4% with a corresponding 
increase in solidity, were shown to be 24%, or twice the pre- 
vious amount. The absolute limit of performance increase by 
this method would be given by the case of negligibly thin 
blades. For this case, the total increase in turning parameter 
would be 50%. 

It was pointed out in this study that increasing the 
performance of cascades by using thinner biades and higher 
solidities would probably reduce the flexibility of performance. 
The effects of operating off-design were not a part of this 
study since the method of generating the cascades produced 
airfoils operating only at the ideal angle of attack. It was 
stated that modifications could be made in the method to con- 


sider the off-design case, but this was not done due to time 


limitations. 
29 







ro 


7 ’ SE a ok 























Peean = 
I i Sa OL ORY me tS <a 
vee ee Oem # ° 
ay ee 
Bee ete tte me ee ome - _ 
210) on ew , ied 
0) Ca ce ee ft ret ot eee oo 





———— | Ns  T monty 
ees eee ee ee ~~. 
ets ime on ‘a 
Ce ee | 

tel) Ved ed See ek od eee 
















—— et ee ed eS ee ee Se Ole 
ete Sl eet eee 
er 6 lle Ot ce eee pee SUITES 
Ce Ot om ee eee Oe be! eee * Sel ot © 
Cee | ete @ lee, oe! ae Cae 
0 ee eee Pe Vk ee vee * eh. 
— eS &  — ——§ Ge ae (ae EA OSES 
- —_ <-> ‘= >= = «€¢ &—————w Sah eas 
ee 











6.0 Suggestions for Future Work 


This investigation represents only a beginning in the 
field of theoretical design. With the basic method and com- 
puter programs available, it should be a simple matter to ex- 
tend the scope of investigations similar to this one. Inves- 
tigations ought to be conducted for a number of average flow 
angles, so as to cover the whole range of possibilities in 
cascade design. Reference 2 for example, has some work on an 
average flow angle of zero, which represents the special case 
of an impulse turbine. There is no reason to prevent con- 
Sidering negative average flow angles either. This would be 
the case of turbines rather than compressors, but the same 
method and computer program would apply. 

One might suppose that since the flow is accelerated 
through a turbine, the boundary layer experiences a negative 
pressure gradient and separation is not a problem. This is 
not exactly true, however. There are actually two pressure 
differences to consider in flow through a cascade. One is the 
effect perpendicular to the stream, which is a result of the 
lift being developed by each blade. The other is the longi- 
tudinal effect parallel to the stream which is the result of 
differences in pressure from the front to the rear of the 
blade row. Flow separation depends on the combination of these 
two effects. Thus, in decelerating cascades, the longitudinal 
effect reduces the allowable lift from each blade. In tur- 


bines, the longitudinal effect increases the allowable lift. 
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If the lift is increased enough to overcome the aid of the 
longitudinal effect, flow separation will still occur. A 
series of studies of this matter could establish the exact 
relation between these two gradients and the occurrence of 
separation. 

The analytical methods of this report would be a help in 
studying the occurrence of separation through solution of the 
boundary layer equations. The detailed pressure distributions 
provided by the computer at least open an avenue toward solu- 
tion of the boundary layer equations by direct integration. 

In any future work, it would be convenient to consider 
the velocity at more points near the leading edge of the blades. 
This would allow more accurate determination of the actual 
maximum velocity, and should eliminate some of the scatter 
found in the data for this investigation. The necessary modi- 
fications to the computer program are not difficult, and only 
lack of time prevented their incorporation in this investiga- 
tion. 

Another area of investigation would be a change in the 
basic transformation which would allow a finite angle at the 
trailing edge of the blades. This would make the transforma- 
tions similar to the Karmann-Trefftz transformations then, 
rather than the Joukowski. Eliminating the cusp or radius at 
the trailing edge then, would perhaps yield blade shapes which 
would be more physically useful. As pointed out earlier, some 


of the blade shapes developed in this investigation, while 
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mathematically useful, were physicaily poor because of the 
rounded trailing edge. These improved transformations would 
be particularly useful in the case of symmetrical impulse 
blading, where a finite angle is desired at both the leading 
and the trailing edge. 

As mentioned previously, it would be interesting and per- 
haps useful to consider the off-design case, where the blades 
are operating at some angle of attack other than the ideal. 

In principle, the method would be to generate the airfoils, 
then go back to the original complex potential and change the 
average flow angle. Keeping the same transformations would 
preserve the blade shape, but the circulation would have to 

be adjusted to keep the Kutta condition. The idea seems simple 
enough, but no doubt there will be many difficulties with the 
details. 

Still further work could be done in investigating other 
Shapes for the camber line of the blades. Generally, the 
camber line for the blades developed in this investigation was 
the slightly flattened arc of a circle. The generation of 
other camber lines involves the expression for the flow poten- 
tial in the near circle plane. In this plane, a single doub- 
let generates the flow contour which transforms into one of 
the blades of the cascade. If instead, a series of doublets 
were used; by varying the strength and distribution of these 
doublets, great control could be exercised over the camber 


line and blade shape of the airfoils. The ultimate in this 
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direction would be where the input is a desired pressure 
distribution, with the proper camber and blade shape being 
automatically produced. The resultant family of cascades 
would be as important a milestone to cascade design as the 


NASA 6-series airfoils were to isolated airfoil design. 
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Figure 6. 
Boundary Layer Loading Frarameter vs. solidity 


For Constant Thickness Airfoils 
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Boundary Layer Loading Frarameter vs. Solidity 


For Constant ‘thickness Airfoils 
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Figure 8. 
Boundary Layer Loading Parameter vs. Solidity 
For Constant ‘thickness airfoils 
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APPENDIX & 


MATHEMATICAL ANALYSIS OF ‘thu tipfnubD ve DeSIGN 


this appendix describes the general mathematical 
development of the method of designing airfoils in cascale. 
wome of the fine details, particularly on iterative proceed- 
ures, have been omitted. However, the entire calculation 
proceedure is given in the runinal program in the next apvpen- 
dix. ‘his, along with the explunations given in the table 
of symbols, should allow analysis of the method to whatever 


extent is necessary. 


Basic nelations in the Near Circle rlane 





the complex potential for the flow around the infinite 
row of closed bodies in the near circle plane may be written 
ice LOrm: 
we pray = cP, + (a + iB) By (2) Ala 


With the corresponding complex velocity: 


-if 


dw =v.-iv.=e + (A + iB) £,(2) Alb 
AG x y 2 
where 
H(z) = K, (x,y) + id, (x,y) = % coth 2 A2a 
#,(2) = Ro(x,y) + idy(x,y) = $5 Fy (2) 
= = + eoth 5 | = ~ + esch® 5 A2b 


By expanding these functions and separating tnem into 
their real and imaginary parts, it is found that: 
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- ~ /Sinh x | A3a 
Hy (xy) ie bes xX - cos y 

Se eee 
Seen) z ( eene xX - cos y ee 


(cosn X 

ae cos = 

Holy) = - 2 (cosh x - cos y Ade 
1 (Sinh x sin 

eae ae: xX - cos y)A aoe 


the expressions for velocity potential, stream function, 


I 


Jo(Xsy) 


ana velocity components now reduce to: 


# 


x cos# +ysinfB + AK, (x,y) = Bd, (x,y) A4a 


== 
Tl 


- x sinB + y cos@ + AJ, (x,y) + BR, (x,y) A4b 
Vo= cos(3 + Alt, (x,y) - BJ, (x,y) A4c 


ca suas) ue Ad, (x,y) - BR, (x,y) A4d 


these relations represent the flow without circulaticn 
about an infinite series of nearly circular bodies s;aced 
at equal intervals S = 277 along the y axis. ‘the contour of 


a typical body is shown in Figure 5. 


Location of Stagnation Points 


bue to the interference caused by th. presence of the 
infinite stack of bodies, the stagnation points on any body 
are displaced slightly from the position they would have 
occupied if the body were isolated in the flow. ror this 
reason, sume computation is necessary to locate these points. 
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vesignating the front and rear stukbnation points by the 
Subscripts 4 and b respectively, anu noting the type of polar 
Symmetry involved, the following fucts are observed. At the 
two Stagn.ution points, both velocity components vanish by 
acefinitivn. also, by reason of the odd symmetry involved, 


<r 


Ma > 7 Xp Vg =~ Vp 
the points a and B lie on the same streamline, ¥, = r V3 : 


and Mee =-~, . However, since 


the only way in which these two conditions on Y can be satis- 
fica aS if Ys —_ Y= QO . since the stream functicn vanish- 
es at the stagnation points, it also vanisnes at all other 
points on the body, since tney are on the sume streamline. 
The velocity pcote:.tial P on the other hand varies from 
a maximum negative value at A to a maxjmum positive value at 
B. the value ¥, then defines the amplitude of the f func- 
tivn on the contour of the body. this amplitude may be 
arbitrarily prescribed, and will then fix the size of the 
near circle contour in relation to the fixed spacing + = 27. 
Let tH) (X Vp) = 44, and so forth. Then; at the pouns 


B from the foregoing discussion: 


oD | re ee We aan. = / 

(, = x, cos+ y, sinf + Ax), - bu), = arbitrary — 

n= -Xp sinf? an gs) cos? + AJ,» + Bi = (U) ADdb 

Vi = cos + AK, - BJn, = 0 Adc 
= ~ ~ Bi = Ada 

Vy = sin (3 AdJ5p ~ Bity, O 2 
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the last two of the above equaxions may be solved simul- 


taneously for the unknown constents a and B with the results: 


= Ry cos 3 + Jon sin 7 
A=  -? —a A6éa 


2 
hop + Jou 
oe sinf IS. cos 3 
2B 2b 


substituting these results into the expressions for f, and ¥ 
yields two equations in wnich the two unknowns are the 
coordinates 72 and Y,; of tne aft stagnation point. The 
equations are trancendental, however, so that an iteration 
technique is necessary for solution. Care was necessary in 
tnis area as all the functions are multi-valued in y, and the 
result could converge to a point diffeiing by some integral 
multiple of 2 in the y dimension from the intended point. 
‘he general method involves starting with low values of 

in effect close to the isolated case wiere the location of 
the stagnation points :s known, una working up to the desired 


value by modest increments. 


igeation of Points on the Contour 


vwnee the location of tue uft stagnation point has been 
found, the constants A and p can be evaluated from equations 
Hoe and A6éb . uhe other points on the comtourscan then be 
located from equations ada and ADb by setting #= 0 and 
asGigning, a suitabl< value to f? >. indisS Computer. progien 
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divides the interval from =, Los, " 


this yields 20 points along each surface of the airfoil, or 


into 20 equal parts. 


40 points in all. As pointed out in the text, it would be 
helpful to subdivide some of these intervals near the leading 
edge to get a better check on the maximum velocity in this 
Crmiuical region. 

ihe equations for the location of these contour points 
are again trancendental, and a similar iterative technique 
is necessary for solution. ounce the coorainates of a point 
are known, the velocity coOtwponents are found from equations 
Adc and abd. dhe final velocity v and tne slope ™ of the 


contour at that point are found from the relations: 


‘f= V = + V 2 A8a 
x y 
V 

ol = arctan a A8b 
x 


At the stagnation points, this relation becomes indeterminate 


and L'Hopital's rule is used for a solution. 


Basic nelations in the Circle slane 


It is well known that the flow around an infinite row 
Morbi trary but identical contours Can aay pmanerpilic we 
mapped into a corresponding flow about a single unit circle. 
ifethe flow is of zero circulation, ang tae civecuicisot eee 
fopecrom the contours is ( then the corresponding flow in 
the circle plene is expressed by the following complex po- 


tential and velocity: 
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4 —". 
we Privy ei ing ~ 1 oth yy at —. 
dw _ : — inane, inf fare = 2 if? /__-2e7 
de = (v,- iv,) + = (2 e aH) + < ote (223) A9b 


uplitting these into real and imaginary parts, and placing r 


= l leads to the result that along the unit circle, 


Me 


A 


0 Al0Oa 
AlOb 


constunt 


1nis verifies that the circle is i: fact a streamline. 


xurther, along tiuis circle: 


f 
V; A f s2fowbe sing + sin sinh & cos & Alld 


cosh 2y" — cos 2 @ 


(coshf’+ cos A 
cosf3 in OOS ae Oe + sin 6° 2 aretan ( Sooo sin © Alla 


cosh¢"- cos@ sinh 


Pancularities anc starvnution points in the Cirele bane 








to find the stagnation point 3B, set Ve =O . ‘then 


equation Allb gives: 
tan @, = tan tanh ¥ ake 


Substituting this into equation alla produces a result which 


can be reduced to the form: 
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: , iemeose \e.. sinl ) | 
qy= cos ln (7 = cosh ) + sinf.2 aretan ( 7 Alda 


wnere 


tT = ) $( cosh 2f\ + cos 2/4) Al3b 


these equations then must be solved for T then %“ for 





assigned values of fs and ¥, which must be the same as the 
values for [s and F previously employed in the near circle 
plane. Since the functions are trancendental, an iterative 


tuchnique is again required. Newtou's metnod is employed, 
q qd Pp 





which involves calculation of the derivative a - bDe- 


? 
noting this by the symbol “a then: 


-81° 


BT ( oq° - cos of)? - 1 


Applying the equations in this form however, requires an 
excessive number of significant figures in the evaluation of 
t . so circumvent this problem, the following subltitution 


was wsed to eliminate i. 


f= cos & +AT AlL5 


the previous equation then become 


2 cos  _ +A , sin? 
P., = cos & inf-2cos +42 I, sinfi 2 arctan cosf “rel 


Al6a 
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the problem of finding ‘ for assigned values of fs and f, 


by wewton's method has no difficulty with regard to signifi- 





cant figures. ‘whe parameter @ and its hyperbolic functions 


may now be found fron: 


sinh ¢“ eee. 2 cos/s +AT) = E ( cosh oy), 


ithe initial estinate for iiewton's method is obtained by 
neglecting JT in the last term of equation Al6a . ‘the 


resulting equution may be solved explicitly for AT. 


A 7 Z cos (S 


oe 
p- (Teeny nen 


Al8a 


where 


cos /9 


these expressions provide an excellent first estimate for 

fo, and in fact for® the special €ase of B= QO , they are 
exact. 

with fs and “a known, the location of any point Oo. 

having an assigned value of P, is defined implicwmily oy 
Cluoaviwn alla. 20 establish the locauion Of poinvs cornes— 
ponding to those in the near circle plane, all that is neces— 
Sacy is to have ti respective values of ?. be the s@me in 
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the two planes. As before, an iterative solution is re- 


quire 


Circulation 





the complex potential and velocity for a purely cir- 
culatory flow of unit circulation in the circle plane are 


as follows: 


if 2 ~ cued 
Aw = - isos m (t= =] al9a 


d (Aw) - 2i cos sinh 1S al 
= 9b 
a3 5 — 2 cosh /*: Suge 


a = re Al9c 


veparating these into real and imaginary parts, it is again 
found that the unit circle r = 1is a streamline, along 
which the velocity distribution is given by 

= sinh 2¢ 420 

Ne eee 
@~ cosh 27° ~ cos 26 

this circulatory component must be evaluated at the sane 
points e. previously establishedsround the periphery of the 
circle. by superposition, the total velocity at these points 


is then given by: 
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Jtagnaticn roints with circulation in Uircle tlane 





Just as the stagnation points a and B are defined by 
the condition v = 0 , so are the stagnation points L and @ 
defined b, the corresponding condition (= 0 . dntroducing 
this condition into the foregoing equations yields the 


essential result 


B : 
@ -.6.2:90.-6 = ja aresin| 72. sinh 2 ¢ cos ©, 
cosh x cos f@ 


A22 


uvnee the points i and It have been located in the circle 
plane, the corresponding values of y. and Ms may be found 
from equation Alla. dintrocucing these values into the near 
circle plane leads to the location of tiie points L and & in 


that prlane by solving equations Ada and A4b. 


wiavoing vetween Circle and Near Vircle lanes _ 





Mapping between the circle cr § plane and the near 
circle or Z2 plane, aS implied in the preceeding discussion, 
is accomplished despite the fact that no explicit function 
hor Zin terms of § or vice-versa is known. +he mapping is 
dene implicitly through tne fact that the non-circulatory 
flows in the two planes must correspond. vr in other words, 
througn the condition that 


w( z) =w (5) A23a 
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bDifferentiating this wives 


dv 
aw} = al. A23b 
aa | az 
la; | 
But 
| = 25¢C 
dz 
and 
adwj= V A23a 
las|° 9 
Hence the mapping derivative becomes 
az = V@ A25e 
[a5 V 


she tens with circulation has been completely defined 
iu the foregoing development for the circle plane. i1n the 
near circle plane, however, only the non-circulatory flow 
has been explicitly defined. «as explained in the text, add- 
ing circulation terms to the complex potential in the near 
circle plane produces unwanted distortions. ‘Therefore, the 
method used is to map the known flow with circulation in the 
circle plane back into the near circle plane using the map- 
ping relations previously established. ‘these mapping rela- 
tions are of course independent of circulation. ‘hus, for 


the final flow with circulation 


We (z) = W ($) =w(S) + Bo Aw (5) A24a 
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Differentiating this gives 





BS a (Aw) 
Moj= (ano was A24d 
az az 
la3 | 


wubstituting in terms of velocities, and introducing proper 
algebraic signs gives finally for the velocity in the near 


circie plane: 


we 


VU = Sess Ale =v + B Ale) 425 


<i 
So 


i 


Basic Kelations in the Cascade Plane 


dhe basic idea of this entire method is to map the near 
circles into a cascade of airfoils in much the ome manner 
that the Joukowski transformation may be used to map a single 
circle into an airfoil. ‘the equations for the transformation 
are identical in form to those already developed for the 
flow in the near circle plane. snowever, all coordinates are 
evaluated with respect to axes shifted slightly from the 
original position. wvenoting these shifted coordinates by 


primes, the mapping function employed is specifically: 


é 


G* = 3 + av seth 4 (A' + iB') F,(z') A25a 
a an _ 3 = -if' 4 ,1 I ' ‘ 
we7 U,- tu, =e + (A' + iB') FA(z') A25b 
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where 


Bi (z') = ¢ coth om A26a 


if 

f 
ai 
Q 

69) 
C 
ro 

NO 


f q — 
ghz") 2 A26b 
Note that the above functions are exactly of tie same form 


aS equations Al and a2. All subsequent results may readily 


be obtained through this analogy. 





Determination of Airfoil vhape 
in order to transform thc near circle into the required 
airfoil shupe, it is necessary to locate the singular points 
rand ,~ of the mapping transformation in the proper relation 
to the stagnation points u and t of the flow. Normally, v 
as located a short distance inside the contour, just oppoSite 
L aa the surface. dimilarly, point ~ lies inside tne con- 
tour on a line normal to the contour at T. xefer to rigure 
5. she significance of tiese locations may be qualitatively 
understood from the fact that as the relative distance i, 
= is made small, the radius of curvature at tne nose of 
the airfoil decreases rapidly, and i: the limit, as € = 0 
tne nose becomes cusped. Similarly, the relutive distance 
sit = Ek controls the tail curvature. in the limit as 
aR =) the tail becomes cusped in tne manner or a Jvoukow- 
Ski airfoil. ‘the origin of coordinates ror the transforma- 


tion is shifted from the original location, the midpoint 














between points A and B, to tne new location, the midpoint 


between points P and QQ. 


Velocity and slope “long the Airfoil 


‘the complex potential including circulation i:. the cas- 


cade and near circle planes may be equated 








x pe f 
wy (§*) = wil) A2Ta 
vifferentiating 
dw 
dw 
Cc = POA A2Tb 
ag (G2 —) (<=) 
dz' az 
T aa Ay 
Let (Se ion a conjugate of velocity A2T¢ 
ag * on airfoil 


(“e) = if a-o A27Ta 
conjugate of velocity 
on meax Camere 


( $3] nes A2Te 
(ge) = ih Peo 


From which, aiter separating amplitudes und phase angles: 





) — : v 

V = 2 = v+B OV, 7 A28a 
tes U 

ie X-0- A28b 
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ArrbNula b 


fortran Program "Cascade" and Uperating instructions 


Jn the following pages is the FuxTiAn program "Cascade" 

used in deriving the data used in this study. input data is 
y means of data cards inserted at the end of the program. 
any number of data cards may be used. the first card indica- 
tes how imany additional data cards are to be rend. 

Lach additional card lists the number of values of the 
parameters , ; Bo ' and K to be computed, their 
initial values, and the increment by which each value increa- 
— Also listed are the print control quantities hlUrn and 
NEAPR, which control the extent of the print out as indica- 
ted in the table of symbols. she computez then calculates 
results for every possible combination of the five parameters, 
prints the results, and then proceeds to the next card. 
the stundard print out contains the five inpv.t parameters, 
the solidity and staggeer angle of the cascade, ashe airfoil 
thickness, and the coordinates, slope, velocity, and boundary 


layer parameter for each of the 42 points of ihe airfoil. 


Condensed #rint vVut 

when surveying large numbers of cascades, as was done 
for this study, it is convenient to condense the data to 
reduce the volume of output. age 85 shows the modifications 
mede to the besic w.rogram to accomplish this. Briefly, att 


ey 


as a 


8 i b ‘ 
= 


= o 2S 
: : 


——— | 


al 


Ab oe som = La ace 

> a a © —=—— —ée ow = 

=< oa. <a et Ga Gey Alem oe 
—- — <a Cosa Gaels , 


= & —_* © =e et C2 eee oe line 


el ed bee - Pr pie 
a ee, 7 —n ——_. —_— mdb 
i «4a a ae “1. o@ —_—— 
"ale Gdn oe” 3 om - —_— -\ «, 
- Be mn me Va ee met hai haul? * gil? 


—aecisee «& Jie ae care i —e 
-enW1 = - & Sees 





—_—= ~*~ a ed a 


=e fe | ol .° =e Sa F == Cm eae 


Ati Ah oii « 


4 = a 


et am ace ‘sel eS ee 


5. ® ° — | j— U 





te _ 





7 - | ‘es 7+ 
—— = ae yw — a mae 
. = =i=e - —_ 
ee es = aS i tha 


— 


of the standard print commends were removed. uhen commands 
to print the heading only were inserted early in the program 
right after changes in sigma. ‘the modified print out then 
ccntains no airfoil data at all except for the maximum value 
of the boundary layer pirameter. vUperation of the program 


is identical with the basic program. 


iterative scheme for Vonstant isahickness Airfoils 





Page 84 shows further modifications made to the basic 
Pte Pa: to produce cascades with the desired airfoil tnick- 
ness. in effect this eliminates the parumeter & since it 
is varied as necessury to produce the required thicknesses. 
ihe initiel value of & is still read from the data card, but 
the column which indicates the number of & to be calculated 
no longer has significance. whe program as written will 
demand thicknesses from 4% to 24% in 4% increments. uther 
values for thickness can be obtained by modifying the appro- 
priate commands. No means of doing this by data cards is 


provided. ‘The print out is of the condensed form described 


previcusly. 
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